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ABSTRACT 

We propose a simple analytical method for estimating the central volume density of prestellar 
molecular cloud cores from their column density profiles. Prestellar cores feature a flat central part of 
the column density and volume density profiles 0/ the same size indicating the existence of a uniform 
density inner region. The size of this region is set by the thermal pressure force which depends only 
on the central volume density and temperature of the core, and can provide a direct measurement 
of the central volume density. Thus a simple length measurement can immediately yield a central 
density estimate independent of any dynamical model for the core and without the need for fitting. 
Using the radius at which the column density is 90% of the central value as an estimate of the size 
of the flat inner part of the column density profile yields an estimate of the central volume density 
within a factor of 2 for well resolved cores. 

Subject headings: ISM: clouds -Methods: analytical - ISM: magnetic fields - ISM: structure - stars: 
formation 



1. INTRODUCTION 

Prestellar molecular cloud cores are the very first 
stages of the star formation process, and as such they 
are the subject of intense observational investigations, 
as their properties can reveal important clues about the 
initial conditions of star formation. Column density pro- 
files of prestellar cores encode information about their 
structure, stability, and evolutionary stage. They are 
measured using various methods, including mm contin- 
uum emission (e.g., Ward-Thompson et al. 1999), mid- 
IR absorption (e.g., Bacmann et al. 2000), dust extinc- 
tion and reddening of background stars in near-IR (e.g., 
Alves et al. 2001; Kandori et al. 2005), and flux measure- 
ments in optically thin lines (e.g., Tafalla et al. 2002). 
These profiles typically feature a uniform density inner 
region (Ward-Thompson et al. 1994; Ward- Thompson et 
al. 1999; Bacmann et al. 2000; Shirley et al. 2000; Schnee 
et al. 2010), followed by a power law decline. 

Detailed fits of various equilibrium (e.g., Bonnor-Ebert 
spheres, Bonnor 1956; Ebert 1955) or dynamically evolv- 
ing models (e.g., Crutcher et al. 1994; Ciolek & Basu 
2000) to column density profiles can provide estimates 
for a variety of core properties. However these detailed 
fits are often computationally expensive to produce, as 
profiles for various parameter values are typically only 
available as numerical solutions. Although analytic pro- 
file fits can alleviate this problem (Dapp & Basu 2009), 
the results of full profile fits are necessarily model and 
geometry dependent. 

Here we propose a simple and model-insensitive way 
to estimate the central value of the volume density of a 
starless molecular cloud core from its column density pro- 
file. Although this technique provides the value of only a 
single quantity, it has the advantage that it is computa- 
tionally inexpensive, as it requires no fitting but only an 
estimate of the size of the uniform density inner region; it 
is insensitive to geometry or amount of magnetic support; 
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and the quantity it estimates is a very important one, as 
the central volume density sets the evolutionary stage 
of the core within the context of any specific dynamical 
model, and so it is necessary for further core modeling, 
including velocity, magnetic field, and chemistry. 

The reason for the robustness of this technique is the 
simplicity of the physics behind the formation of the uni- 
form density inner region: the fiat inner region of the col- 
umn density profile corresponds to the scale within which 
thermal pressure erases any inhomogeneities, the thermal 
lengthscale. It is therefore expected to be present inde- 
pendently of the details of the core dynamics, such as 
the amount of magnetic support of the core. Its extent 
is also not very sensitive to geometry effects or core ori- 
entation, as in the uniform density inner region pressure 
forces, even in magnetic cores, are much more prominent 
than further out, and the geometry is closer to spherical 
than at larger scales. Since the thermal lengthscale de- 
pends only on the value of the central volume density of 
the core and the core temperature, which, in prestellar 
cores, typically varies by only a factor of ^ 2, a measure- 
ment of the thermal lengthscale, using the size of the flat 
inner region of the column density profile, can yield an 
estimate of the central volume density with comparable 
(factor of ^ 2) accuracy. With independent estimates of 
the kinetic temperature (e.g. from line intensity ratios), 
the accuracy on the central volume density can improve 
further. 

In [2] we describe our proposed technique in detail, we 
demonstrate its robustness, and we discuss systematic 
uncertainties entering estimates of the central volume 
density obtained this way. In ^ we apply our method 
and obtain estimates for the central volume density of 
the starless core B68 and the prestellar core L1544. We 
discuss our findings in f|4l 

2. METHOD 

We first give a simple demonstration of the inner flat- 
ness of the column density and volume density profiles, 
and the similarity of the observed extent of the uniform- 
density inner region in the volume density and column 
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Fig. 1. — Radial profiles of the number volume density n (panel a); and column density (panel b). Dotted lines: non-magnetic model; 
solid lines: magnetic model. Each curve corresponds to a snapshot in time when the volume density is an order of magnitude higher than 
the previous one. Time elapsed per snapshot in magnetic model: 4.85 Myr, 5.61Myr, 5.74Myr, 5.77My; in nonmagnetic model: 0.89 Myr, 
1.02Myr, 1.05Myr, l.OGMyr. The fifth (lowest-density) snapshot in both cases is the initial configuration of each model. The thermal 
lengthscale (Eq. [1} is indicated with a circle for the magnetic model and a square for the nonmagnetic model. 



density profiles, for two dynamical models of prestellar 
molecular cloud cores. The first is a spherically sym- 
metric non-magnetic dynamical model. The second is 
an axially symmetric magnetic model, where the initial 
mass-to-flux ratio is equal to the critical value for col- 
lapse (Mouschovias & Spitzer 1976). Details for these 
simulations are given in Tassis et al. (2010, in prepara- 
tion); the two models discussed here are the "reference" 
non-magnetic and magnetic models in that work, respec- 
tively. 

In Fig. [1] we plot radial profiles of the volume density 
(left panel) and the column density (right panel), for the 
non-magnetic (dotted lines) and magnetic (solid lines) 
models, for different time snapshots. Each snapshot cor- 
responds to a time when the central volume density is an 
order of magnitude higher than the previous one. The 
time elapsed between the initialization of each simulation 
and each snapshot are given in the caption. 

The radial profiles of the volume density and the col- 
umn density have similar properties, and feature a flat 
inner region, the boundary of which is indicated with a 
circle for the magnetic model and a square for the non- 
magnetic model, followed by a power-law decline with 
increasing radius. The angular size of the column den- 
sity profile flat inner region is observationally measur- 
able through a variety of techniques (see discussion in 
iJT]). Combined with an estimate for the distance of the 
core, such observations can yield the physical size of the 
column density profile flat inner region, which, however, 
is of the same size as the volume density profile flat inner 
region, as can be directly seen by comparing the left and 
the right panels of Fig. [TJ The latter corresponds to the 
thermal lengthscale (Mouschovias 1991), which is equiv- 
alent to the critical Bonnor-Ebert sphere (Bonnor 1956; 
Ebert 1955), and it depends on the temperature and the 
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Fig. 2. — Central density as a function of thermal lengthscale 
(Mouschovias 1991) for four different core temperatures. Solid line: 
T = lOK; dotted line: 8K; dashed line: 15K; dot-dashed line: 
20K. Because observed core temperatures only span a small range, 
observing Xt can give us the central volume density. 

central density of the core: 

/ y X 1/2 /io3cm-3V/2 ^ 

^--«^io(iok) (^^) (1) 

Because the temperature does not vary much between 
prestellar cores, measuring At is a proxy to measuring 
the central volume density of a collapsing core. The very 
modest (factor of 2) variation of Uc for a given value 
of At due to the corresponding factor of ^ 2 spread in 
possible core temperatures is shown in Fig. [21 If the core 
temperature can also be independently measured, then 
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Fig. 3. — Fractional error, for various dynamical models and 
beam sizes, in the calculation of the size of the thermal length- 
scale from the column density profile, when using the criterion 
N{\t) = 0.9N(0). Thin dotted line; non-magnetic spherical 
model; thin solid lines: magnetic models, viewed edge-on, symbols 
indicate the initial mass to magnetic flux ratio (diamond: subcrit- 
ical, magnetically supported; x: critical; -|-: supercritical). Dashed 
curves are magnetic models viewed face-on (symbols as for thin 
solid curves). Normal solid and dotted lines correspond to a beam 
size half the angular size of \t, for magnetic and non-magnetic 
models respectively. Thick solid and dotted lines lines correspond 
to a beam size equal to At , for magnetic and non-magnetic models 
respectively. 



the accuracy with which Uc can be obtained through a 
measurement of improves even further. Equation ([T]) 
aUows one to obtain an actual volume density measure- 
ment directly from a size measurement. This result is 
independent of dynamical model and core geometry, as 
it is a very basic consequence of thermal pressure. 

In practice, we recommend determining At from the 
column density profile N{r) through the criterion 



N{Xt) = 0.9N{0) . 



(2) 



This criterion is independent of the absolute calibration 
of the column density profile and is not affected by any 
uncertainties in this calibration. The fractional value of 
0.9 at the "edge" of the uniform-density inner region is 
optimized so as to minimize sensitivity to the details of 
the dynamical model, beam size, and core orientation 
with respect to the line of sight. 

The robustness of this criterion, as well as a quan- 
titative estimate of the systematic uncertainties in the 
determination of Xt and entering through these con- 
siderations, are demonstrated in Figs. |3] and 01 Figure 
[3] shows the fractional error in the determination of At 
when the criterion of Eq. ^ is used, for various dynam- 
ical models and beam sizes, detailed in the caption. In 
magnetic cores the criterion of Eq. ([2]) tends to over- 
estimate At at lower central densities, while the trend 
is reversed at higher central densities (and later times). 
The physical reason behind this systematic trend is the 
extra magnetic support perpendicular to the field lines 
in magnetic cores, which is decreased at higher densities. 
The beam size does not significantly degrade the accu- 
racy of the At determination, as long as the flat part of 
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Fig. 4. — Estimated vs true value of the central volume density. 
Line types and symbols as in Fig. |3] The long dotted line denotes 



the profile is resolved. Even with a beam size compara- 
ble to the size of At (thick lines in Fig. [3]) , the fractional 
error does not exceed 50%. The difference between face- 
on (dashed lines) and edge-on (solid lines) orientations 
of magnetic cores is small. 

The corresponding error introduced to the estimate of 
He through these systematic uncertainties affecting the 
determination of At is shown in Fig. |4l If we do not con- 
sider marginally resolved inner core regions (i.e. beam 
sizes ~ At; thick lines), the deviation from nc,truc is al- 
ways smaller than a factor of 2. Even with the inclusion 
of the large beam size cases, the deviation does not ex- 
ceed a factor of 3. 

3. APPLICATION 

As an example of our proposed recipe, we apply the 
simple technique described above in two well-studied 
molecular cloud cores, B68 and L1544. 

For L1544 we can read off of both the 450 and the 850 fi 
m profiles of Shirley et al. (2000) through a simple linear 
interpolation between points that the column density (in- 
tensity) declines to 90% its central value at a distance of 
rgo ~ 1400 AU from the center of the core. Using a tem- 
perature of 12 K (e.g., Williams et al. 1999), we obtain 
an estimate of the central volume density of nc ^ 2 x 10^ 
cm~^. This is significantly higher than the estimate of 
Williams et al. (1999), however it is consistent with the 
higher estimate of ~ 10^ cm~^ of Ward- Thompson et al. 
(1999) and it is in excellent agreement with the value ob- 
tained by Ciolek & Basu (2000) through detailed fitting 
of a magnetic model to this core. 

For B68 we use the data of Alves et al. (2001) to obtain 
an estimate of rgo ~ 2700 AU. For T w lOK (Hatzel et 
al. 2002; Lai et al. 2003) we obtain ~ 5 x 10^ cui^^. 
This value is consistent, within systematic uncertainties 
discussed above, with the value of 3 x 10^ obtained by 
Dapp & Basu (2009) by detailed fitting to their analytical 
approximation of a Bonnor-Ebert profile. 

4. DISCUSSION 
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We have shown that, through a very simple analytic 
recipe that does not require fitting, it is possible to obtain 
an estimate of the central volume density of a starless 
molecular cloud core from its column density profile. 

For the optimal application of this technique, the beam 
size used should be small enough so that the flat part 
of the column density profile is well resolved, and the 
temperature of the core should be independently known. 
Under these conditions, through a simple measurement 
of the extent of the uniform column density region of 
the profile, we can estimate the value of the thermal 
lengthscale and from there the central volume density, a 
valuable input for further modeling, and an evolutionary 
"clock" for any specific dynamical model. This proce- 
dure returns a model-insensitive estimate for the central 
volume density within a factor of two. 

Uncertainties in the estimate of the central density in 
addition to variations in the shape of the profile due 
to the varying amount of magnetic support can arise 
through the following effects. 

(i) If the temperature of the core under consideration 
is not independently known, additional systematic 

uncertainty, proportional to the uncertainty in the 
value of T, will enter our final estimate of Hc- 

(ii) Although the temperature gradient in the inner flat 
part of the profile is expected to be small (e.g., 
Tafalla et al. 2002), any deviations from isother- 
mality will similarly induce uncertainties in the de- 
termination of the density profile and the estimate 

of He- 

(iii) Deviations of cores from spherical symmetry will 
impact the shape of the profile even close to the 
flat inner part and can affect the estimate of ric. 
Cores typically do have aspect ratios of ~ 2, and 
although we have considered here flattened cores in 
the case of magnetic support, it is possible that as- 
pherical cores may arise for other reasons; however 
the uncertainty induced on Uc in such cases is not 
expected to exceed that which we have explicitly 
calculated for cores deviating from sphericity due 
to varying amounts of magnetic support, which also 
exhibit aspect ratios of 2 or higher. 
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(iv) The calculation of a central volume density requires 
the knowledge of a physical, rather than angular, 
size of the flat part of the column density profile, 
and therefore any uncertainties in the distance of 
the cloud will also propagate to the estimate of the 
central volume density. 

Depending on the detailed conditions in the parent 
coud and the details of the formation mechannism, the 
structure of the core at radii significantly larger than the 
central fiat part of the column density profile can differ 
appreciably from the profiles presented here. However, 
the density profile in the most central parts of the core 
that we consider here, and at central densities in excess of 
~ 10^ cm~'^, has lost, to a large extent, memory of the de- 
tails of the core formation mechanism and of the physical 
conditions at the envelope. This has been shown by, e.g., 
Eng (2002) for the presence of MHD waves in the parent 
cloud, Tassis & Mousehovias (2007) for varying amounts 
of initial magnetic support, and Gong & Ostriker (2011) 
in the case of cores created by colliding turbulent flows. 
The latter, once collapse starts, exhibit density profiles 
compatible with models of pure hydrodynamical collapse 
(Larson- Penston solutions). For this reason, and as long 
as the analysis is performed very close to the central part 
of the core, our results are relatively robust to the nature 
of the processes which are operating at larger scales and 
which are responsible for the formation of the core. 

Although more detailed modeling can yield more ac- 
curate estimates of the volume density and an array of 
other parameters, the remarkable simplicity of our pro- 
posed technique and its insensitivity to the details of the 
adopted dynamical model render it a fast and useful tool. 

We have verified that the results obtained in well- 
studied molecular cloud cores LI 544 and B68 are con- 
sistent with other, more detailed estimates available in 
the literature. 
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